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ABSTRACT 

The brightness temperature of Titan at 3 mm wavelength is anmnd 200 K, 
according to Ulich, Conklin, and Dickel (1978). Although an earlier 
measurement by Briggs is much colder, we ack^t 200 K as the surface 
temperature and build an atmospheric model with a surface pressure of 21 
bars. CH 4 clouds form between 100 and 120 km altitude. The visual limb 
is near 200 km. The methane mixing ratio is 0.25% above the clouds and 
7% below; the dominant gas is assumed to he Nj. - 'le thermal opacity is 
due to pressure-induced absorption in Nj and a trace (0.5%) c ** **- with 
some help from cloud particles; unit opacity is reached at 600 mbar, 1 lU km 
from the surface. The radius of the solid body in this model is 2700 km, in 
reasonable agreement with 2600 km obtained if the density is the same as 
that of Ganymede and Callisto. Deeper atmospheres can be obtained if the 
temperature gradient is subadiabatic or greater CH 4 abundance is assumed. 


INTRODUCTION 


It has been clear for some years tiiat a deep, cloudy atmosphere, with a surface 
pressure of many bars, was not ruled out for Titan by any existing data (Lewis and 
Piinn, 1973; Pollack, 1973). On the other hand, no data required it. Recently, 
however, Ulich et dL (1978) have obtained a 3-mm br^htness temperature of 200 K. 

It is, therefore, worthwhile to eiqilore the properties of a deep atm isphere in some 
detati. The model has quite a bit in common with the established picture of Venus, 
with a surface temperature 3-4 times the eCGective temperature, and a dense cloud 
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layer tens of km deep. Because of the lower gravity, however, this depth is achieved 
with a lower surface pressure, only 21 bars. 

Since the next three sections are devoted to a variety of detailed discussions, 
we refer ahead to the final model illustrated in Figure 3. There is a clear layer 
100 km deep composed of N2 ^ ^ Dense methane clouds form at 

100 km and extend up another 10-20 km. Their top is diflfuse enough to have an impor- 
tant effect on the formation of the observed absorption bands. We do not treat the 
question of how' much solar radiation penetrates to the surface to drive the greenhouse, 
but w« can appeal to the anal(^’ with Venus to suggest that the mechanism is at least 
plausible. 

The model presented here is far from unique, but we feel it to be conservative, 
in the sense that any plausible variant with a 200 K surface would have a still deeper 
atmosphere and greate’- surface pressure. An absorbing layer at the 200 K level is (me 
possibilify; another has more methane vapor, which increases die scale height and 
gives a deeper region of small (wet adiabatic) lapse rate. Or the lapse rate could be 
less than the adiabatic. 

A summary of numerical data for Titan is given in the Appendix. Earlier 
review's appear in the books edited by Hunten (1974) and Bums (1977). We make no 
attempt to summarize this material, but some of it is surveyed and updated in the pre- 
ceding article by Caldwell. 


THE RADIO DATA 


The brightness temperature T^ deduced from a radio flux measurement is 
inversely proportional to the sqture of the assumed radius. Table 1 shows die data of 
Ulich et al (see also Conklin el al., 1977) and Briggs (1974) as given and as converted 
to a radius of 2700 km, justified below. 


Table 1. Radio Data 


Source 

Wavelength 

(mm) 

tb 

(K) 

Radius 

(km) 

T„ (2700) 
(K) 

Ulich et al. 

3 


2900 

231 

Briggs 

37 

103 

2530 

115 
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Although each result has fairly lai'ge error bars, tiiey do not overlap. Sin^ we do not 
know Ik>w to resolve the discrepancy, wo take a round value of 200 K for the present 
work, assuming it to be the surface temperature. The bias above the men of the data 
is partly to account for an emL««ivi^’ less than 'mi>.*. radiatii^ layer could also be 
within at3ncs|>h;s^e if thCie is enough mtcrow’avc opacity, dt^ perhaps to NH^. For 
convenience, however, we shall use the term * surface". A lower temperature, if 
required by later data, is readily accommodated by inserting a new "surface" at the 
appropriate height. 

COMPOSITION OF THE VISIBLE ATMOSPHERE 

Quantitative interpretatiim of the "red" bands of CH4 (which extend from 4410 
to 9000 A) has only recently become possible through the laboratory work of Lutz. 

Owen, and Cess (1976) (hereafter. LOC), Fink, Benner, and Dick (1977), and Giver 
(1978). For practical purposes, the absorption at each waveleng^ is purely exponen- 
tial and independent of pressure; an observation of Titan, ttierefore. gives a methane 
abundance, and LOC obtained SO m-A (meter-Amagat). Lines of the band are 
pressure-dependent in a known way; LOC used a result by Trafton (19? 5b). with their 
methane abundance, to obtain a total abundance of 21 km-A (if the major gas is N2); ttie 
corresponding pressure is 300 nibar for a gravity’ of 117 cm and the CH^ mixing 
ratio 0.4^'. There is no direct evidence that the gas is N2. but we feel the indirect 
arguments for it are strong (Hunten, 1973, 1977). 

This analysis assumes that the visible atmosphere is a clear gas above a dis- 
crete cloud top; the reflecting layer model or RLM. If, on the other hand, the cloud 
top is very’ fuzzy, the homogeneous scattering model or HSM is moi'c appropriate. 
Indeed, Trafton (1975a) has already clearly demonstrated that the RLM does not fit the 
data; the analy’sis below’ merely substitutes laboratory data for the Saturn data he was 
forced to uce. A crude HSM analysis has kindly been carried out by’ L. Wallace, as 
illustrated in Figure 1. The residual intensities, as a fracticn of the nearby continuum, 
were estimated for several wavelengths from the spectra published by Trafton (1975a), 
and were plotted against the laboratory absorption coefficients of Giver (1978). It can 
hardly be claimed that the HSM. represented by' the solid curves, is a good fit to the 
points, but at least it is much better than the RLM, w’hich is simply an exponential. 
Some of the discrepancies may derive from the assumed location of the continuum at 
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longer waveienglte: it may be tiiat the continuum is never reached. The best fit lies 
between the two curves of Figure 1; we adopt a specific methane abundance M of 
150 m-A per mean free path. 

To interpret tte 31^3 line, we use flie curve of growth (Oieir Figure 12) derived 
by Wallace and Hunten (1978), which takes account of the varying pressure through the 
scatterii^ atmosphere. With the specific abundance and equivalent width as input, we 
get a value for pj^, the pressure at unit scattering optical depth. For this first rough 
analysis, we divided Traflon’s (197%) equivalent width of 1.0 A, and the laboratciy 
strength, by 2, since the line can be approximated as a doublet with only slight overlap. 
The broadenii^ coefficient at 1 bar is a® = C. 16S cm“^ for N2 at 100 K (Damton and 
Margolis, 1973). ‘Ihe '.esuU is p^ = 900 mbar; the abundance corresponding is 
63 km-A which, like die CK^ specific abundance, is referred to unit mean free patii. 
The ratio, 0.24%, is dierefore, the CH4 mixing ratio according to the HSM. It makes 
little difference whether we use die RLM or die HSM, but we shall adopt the latter 
value, which is clearly uncertain by at least a factor of 2, 
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Methane photolj'sis and escape of the resulting H2 gives a predictable H2 mixing 
ratio of 0.47*7 (Hunten, 1973, 1977). Our picture of the visible atmosphere is 
therefore: 

- A gas, predominantly K2 containii^ 0.25^ CIl^ and 0. 5fx H2: 

- A cloud (argued below to be mainly Ikjuid CH4) uitii a very diffuse top, and 
unit optical depth at 900 rabar; 

- A haze of *'Axel dust" (Danielson et oL, 1978; Caldwell, 1977) extending well 
into the stratosphere. Since it is probably a mixture of methane photolysis products, 
it is likely to dissolve readily in the methane cloud droplets and give them also a dark 
color. 

Elliot, Veverka, and Goguen (1975) have obtained a radius of 2900 km by analysis 
of a lunar occultation. It is hno^^m that a scattering atmosphere closel}' resembles a 
Lambert sphere, ttte model that was used (Harris, 1961). But what pressure level does 
the radius refer to? A ray grazing the limb encounters unit optical depth due to mole> 
cular scattering alone at a pressure of 50 mbar. hi Titan's impure atmosphere, th^ level 
might be higher by 1 or 2 scale heights, or the pressure as low as 10 mbar. About 4 
scale heights, or 100 km, separate ttiis level from the 900 mbar level; the latter is, 
therefore, close to 2?00 km radius, a value we henceforth adopt. 


FAR-INFRARED OPACITIES 

It is now necessary to tie the optical structure to the thermal structure. As on 
I’enus, a high surface temperature can be sustained only 1^ a large opacity, probably 
greater than 100, through the whole thermal infrared. Pollack (1973) has laid much of 
the groundviork for study of a Titan greenhouse. We have simply adopted his illustra- 
tive result for pressure-induced absorption in II2, scaling tte factor 2. 67 (Kiss, 
Gush, and Welsh, 1959) appropriate for collisions with N2 instead of H2. Nitrogen 
opacity (Bosomut>rth and Gush, 1965) dominates at the lower frequencies below 
330 cm~l, and totally suamps the CII4 absorption which has a similar shape. 
Bosomworth and Gush give N2 data only at 298 K; we have omitted the temperature 
correction to both width and strer*]^h of the absorption, discussed by them and also Ity 
Pollack (1973). The adopted curve is shown in Figure 2. 

As found below, the gases alone have plenty of opacity to blanket the warm sur- 
face. But it is known that the structure due to the H2 S(0) and S(l) lines does not appear 
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in Titan’s thermal emission specti-um (I’oIInck, 1973; Ixnv and Rieke, 1974; Hunten, 
1977). It is reasonable to suppose that tlu' eloiid particles provide enough opacit\' above 
350 cm~l tu smooth out the structure. Liquid methane ahmt' would not Ix' adequate, 
but the dissolved hydrocarbons already roterred to might give enough absorption at al! 
fiequencies. Following Hanielson et al. (1973). w have adopted the curve shown in 
Figure 2, with opacih' proportiemal to frequency, I’se of the IISM implies uniform mix- 
ing of particles and gas: if, however, the particle mixing ratio increases downwards, 
the corresponding absoi^ption can be cnidely repivsented in thr' same terms as the N 2 
absorption. 

We now form the Rosseland m.ean of the data in Figuiv 2 (e.g. , Sohwarzschild, 

1958): 


_L f <!5_ 

where x = hr/kT and A is the absoiption coefficient. Without cloud opacit\-, the 
result is 


( 1 ) 


A r 


1.87 X 10"' cm 
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with clouds included, we have 

A = 3. 10 X 10"^ cm'^ A"^. 
c 

The level having temperature the effective temperature, is taken at the depth where 
the optical depth r is unity. For pressure-induc^xl absorption, 

T - ^ wVj). m 

where H is the scale height and w(N2) is the abundance. Solving for the latter, and 
converting to pressure, we find 

p = 790 mbar (no cloud) 

= 610 mbar (with cloud) 

We adopt the latter estimate and give it the equilibrium temperature 77 K. It is slightly 
above the 900 mb level that has unit cloud o{^cal depth in the red. For tlie strato- 
sphere we simply accept the thermal-inversion model of Caldwell (1977 and preceding 
chapter of this volume). 

THE DEEPER ATMOSPHERE 

We shall, somewhat loosely, use the name "tropopause" for the point (610 mbar, 
77 K> just derived. For the deeper structure we assume an adiabat, making allowance 
for tte latent heat of condensation of the methane that forms the cloud. 

Sagan (1969) has given a conveniently parameterized discussion of radiative and 
convective temperature gradients, and finds them so nearly the same that the tropo- 
pause location is often very sensitive to small changes. Conversely, the profile is 
insensitive to just where the tropopause occurs. 

As a first approximation we neglect latent heat and find the pressure that cor- 
responds to 200 K. The adiabatic lapse rate for N_ is F = pg./C = 1. 16 K/km and the 
scale-height gradient is ^ = -1/3.5. With the standard expression (e.g. , Hunten, 

1971) for the larometric law, 

we can substitute ttie ratio of temperatures (200/77) for H/H^ and obtain a pressure 
ratio of 28. 2. This estimate of the surface pressure is tiierefore, 0. 61 x 28. 2 = 17. 2 
bars, .'\fber allowance for condensation as described below, a better estimate is 19.7 
bars of N2. A possible surface material is methane clathrate hydrate, whose vapor 
pressure is 1.38 bars at 200 K (see Appendix). The atmospheric composition at the 
surface would thus be about 6. 6% CH4, 93^ N2, which gives a total pressure 
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Pg = 21. 1 bars, as well as plenty of material to condense into the clouds. Since the 
vapor pressure of CH 4 is about 2 orders of magnibde greater than tiiat of clathrate, 
there will be an exteikled clear region above the surface. Once ti» cloud is reached, 
the dry adiabat used for the preliminary estimate is not satisfactory. Various expres- 
sions for the wet lapse rate appear in the literature (Brunt, 1933, 1939; Danielson 
et a/., 1977), but the best seems to be that of I^ker (1963). In practice the differences 
are minor, especially since the fate of the condensate is uncertain. If it falls out or 
remains at the height where it condenses, it? specific heat should be omitted, and the 
curve is called a "pseudoadiabat". If the ccmdensate moves with the gas. the situaticm 
is truly adiabatic. Lasker's expressicm for Qie pseudo-adiabat is 


a _ d In T _ 1 Y(1 + ot) 

^ ' d In p “ (Cp/R) + y a(l + a) " 


(4) 


Ifere ^ is the same quantity used in (3), vdiich, luiwever, is valid only over a r^on 
Aj^herc p is constant. Hie specific beat pertaim to the mixture of gas and vapor; 
y is the mixing ratio fCH 4 ] / [N 2 ] ; and -a is ffie exponent in the vapor-pressure rela- 
tion. For CH 4 , a= 1024/T (see Ai^miix). As suggested by Lasker, a suitable method 
is to find the (T,p) relation frmn (4) and then obtain tiie height scale from the hydro- 
static equation. Where is indepeiKlent of height, (3) may be used, and it gives useful 
accuracy even where p is slowly varyii^. In ffie present case, y varies from 0.070 to 
essentially zero and a ^ 12. 5; - jS is 0. 282 below the cloud, and drops to 0. 127 at the 
cloud base, back to 0.286 above the cloud. 

Figure 3 shows the model. The vertical scale is linear in tiie logarithm of 
pressure, with an auxiliary Nonlinear) height scale. The Danielson-Caldwell strato- 
spheric model does not specify a pressure scale; the one adopted is by analogy wiffi ffie 
Jupiter model of Orton (1977). 

The cloud base is at 103 km and has a methane partial pressure of 70 mbar at a 
total pressure of 1080 mbar. The potential cloud mass is the amount contained in a scale 
height, or 600 g cm~^. In the form of 10 pm diameter droplets, this amount would have 
a scattering optical thickness of over 10 ^. Thus, most of it could precipitate out and 
still leave a very respectable cloud. The temperature is 87 K, somewhat below the 
freezing point (89 K); the form may ther^ore be either crystals or supercooled drop- 
lets. The level found in Section 3 for unit 0 |Aical depth in the "red" is a few km above 
the base, a reasonable location. 
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OISCITSSION 


As pointed out in the Intrt^uotion, the pi'esent model is conservative in the 
sense that any other that tcives a "surface** tem^ierature of 200 K roquiix»s an even 
deei>er atmospheix'. A totally different aiiiument is to estimate the radius of solid 
Titan from its known mass and a reasonable me:m densih'. For the latter, ne take 
Callisto ;uid Cuxnymedo as ahaKi»is. weijrhtins the densities of Morrison cf ai. (1977) 
inversely by the stated errors. The i\?sult (1. S9 g cm~3) combines with a mass of 
1.40 X 10--'’ jj to jjive a radius of 'JOOo km. A 10\ chaiifre of densih' reflects into a 
radius change of km. Titan is likely to be less dense, if anything, than any Galilean 
satellite; if so, a reasonable range for the radius is 2000-2700 km, which just includes 
the value from tl»e atmospheric model. 

The vapor pressure of Nlfg dissolved in ice is roughly I inbar at 200 K. Tne 
mixing ratio could thus be around 5 x t0~'*. The optica! depth of this amount is about 
S at 3 mm wavelength, according to the approximate treatment of Field (1959). Per- 
haps the ammonia is kept out of the atmospliere by a layer of hydrocarbons on the sur- 
face. 

Much of this paper is built on the assumption that the surface temperatui'e is 
200 K, as suggested by the microw:ivo obseiwations. The issue of ammonia opacity 
only serves to point up the crucial nature of such data, and the importance of data at 
other wa\x* lengths. The dis;igreement shown in Table 1 is particularly disturbing, and 
a resolution is badly needed. The measurements are, however, so difficult that 
progivss is sure to Ix' slow. Ammonia opacity cannot be called on, because it is nearly 
the same at the two wawlengths. The preliminary results reported by Caldwell 
at the end of the preceding chapU'r support a lowv'r temperature, although they would 
still jx'rmit it to tx’ alxne 150 K. Howwer, a temperature of, say, 100 K c:m be 
accommodated b\ placitig tiie surface at 2 ixirs in' v“<5 km in Figure 1. rhc model is 
unchanged at greater heights. 


Added after the Workshop: 

Shortly after the Workshop, a preprint by Podolak and Giver (197 S) became 
available. Ihis work, an extension of Pixlolak and Oanielson (1977), adopts as its model 
an atmosphere of pure CII4. 2 km- A deep, with a layer of Axel dust (x*cupying the top 
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50 m-A. This dust is nearly opaque in the blue, grading rather suddenly in the red to 
essentially transparent in the infrared. The surface (or cloud top) has an albedo of 
50-60‘I. An excellent fit is obtained to the observed CH 4 bands from 5000 to 10,000 A. 
Although the 81/3 band is not discussed, the model is consistent with Trafton's analysis 
for a reflecting-layer model, which requires 1. 5 km-A if the atmosphere is pure CH 4 . 
At this and longer wavelengths, the Podolak-Giver model reduces to a reflecting layer. 

How might a choice be made between this model and the one described here? 
They were both fitted to the same data; thus we require some additional information. 
Two recent lines of evidence can in fact be brought forward: ( 1 ) the work of Rages and 
Pollack reported in Chapter 11 of this volume; (2) spectra of the 1-3 gmi region. 

Rages and Pollack find the mean particle radius to be 0 . 2-0.4 pm, and strongly 
exclude the value (slightly under 0. 1 pm) required by Giver and Podolak to obtain the 
necessary optical properties. 

The infrared spectra, shown in Figure 4, were obtained U. Fink and 
H. Larson in 1975, and are reproduced here 1^' their kind permission. 

A detailed analysis has not yet been done, but two facts stand out clearly; 

( 1 ) The Titan spectrum is utterly unlike the laboratory spectrum of 1 . 5 km-A 
of methane. 

(2) The closest analc^ to Titan in this spectral r^ion is Neptune. 

The laboratory spectrum is the analog of a clear atmosphere of pure metiiane 
(the Podolak-Giver model at these vrave lengths). Neptune is an analog of a deep scat- 
tering atmosphere with a small fraction of methane. Although the major gas is 
assumed to be H 2 -He on Neptune and N 2 on Titan, the important thing is the scale 
height, not the composition, and the scale heights are rather similar. These argu- 
ments strongly suggest that the deep, N 2 -Cii^, atmosphere is preferred over the 
pure-CH^ one. A detailed analysis is highly desirable, but in this spectral region 
it is very difficult because of the fine structure of the bands. 

A direct comparison of a laboratory spectrum with Titan may be misleading 
because of the gross differences in the temperatures and pressures at which the CH 4 
is obseived. It has been shown by Ramaprasad et al. (1978) in a study of liquid CH 4 
that band shapes are somewhat different at 90 K than they are at a room temperature. 
For example, a band may be seen at 9170 A in the spectra of the cold liquid and of 
' Titan which (s completely masked by the wing of the strong band at 8890 A. The 
sensitivity of the populations of high J levels to temperature must be accounted for 
before a truly meaningful comparison can be made. 
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Figurt 4. 1-2 5 fim spectra obtained tn by U . Fink and H. Larson, reproduced by their kind 

permission. A solar-type star shous the atmospheric transmission. 
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Clearly, all this evidence refers only to the region that can be probed 
optically, and bears only on the composition of tliis legion. The issue of the surface 
temperature and pressure is totally indepei«lent. 

Atreya et al. (1978) have discussed in detail the production of by photolysis 
of NHg, and have found that a surface partial pressure of 14-19 bars could reasonably 
be generated over the age of the solar system. 
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APPENDIX 

Table 2 gives various quantities for three radial distances r corresponding to 
heights of 0, 100, and 200 km in the model. They ai*e the acceleration of gravity g, 
the scale height for pure at a refei'ence temperature of 100 K, tlie base pressure 
for an Ng abundance of 1 km-.\magat, and the dr>- adiabatic lapse rate for Ng. 

Table 2. Quantities for Three Ifadial Distances 


r (km) 

2700 

2800 

2900 

g (cm s"“) 

125 

117 

109 

H (km) 

23.8 

25 . 0 

27.2 

p/N (mb/ km- A) 

15. 6 

14.6 

13.6 

r (dry) (K/km) 

1.20 

1.13 

1.05 


Table 3 gives constants in tlie vapor-pressure equation, p = A exp (-B/T), 
from Delsemme and Wenger (1970) and the H:mdbook of Chemistry and Physics. The 
melting point of methane is 89 K. 

Table 3. Constants in the Vajwr- Pressure Equation 


Form 

A(bar) 

B =oT 

Clathrate 

75,300 

2182 

Liquid 

9,714 

1024 

SoUd 

59,662 

1190 
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